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We analyse nonlinear wave propagation and cascaded self-focusing due to second-harmonie génération 
in Fibbonacci optical superlattices and introduce a novel concept of nonlinear physics, the quasiperiodic 
soliton, which describes spatially localized self-trapping of a quasiperiodic wave. We point out a link 
between the quasiperiodic soliton and partially incohérent spatial solitary waves recently generated 
experimentally. 
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For many years, solitary waves (or solitons) have 
been considered as cohérent localized modes of nonlin- 
ear Systems, with particle-like dynamics quite dissimilar 
to the irregular and stochastic behaviour observed for 
chaotic Systems ||]. Howcver, about 20 years ago Akira 
Hasegawa, while developing a statistical description of 
the dynamics of an ensemble of plane waves in nonlin- 
ear strongly dispersive plasmas, suggested the concept 
of an incohérent temporal soliton, a localised envelope 
of random phase waves ||. Because of the relatively 
high powers required for generating self- localised random 
waves, this notion remained a theoretical curiosity until 
recently, when the possibility to generate spatial optical 
solitons by a partially incohérent source was discovered 
in a photorefractive médium Q| , known to exhibit strong 
nonlinear effects at low powers. 

The concept of incohérent solitons can be compared 
with a différent problcm: the propagation of a soliton 
through a spatially disordered médium. Indeed, due to 
random scattering on defects, the phases of the individual 
components forming a soliton expérience random fluctu- 
ations, and the soliton itself becomes partially incohér- 
ent in space and time. For a low-amplitude wave (linear 
régime) spatial incohérence is known to lead to a fast 
decay. As a resuit, the transmission coefficient vanishes 
cxponentially with the length of the System, the phe- 
nomenon known as Anderson localisation Q. However, 
for large amplitudes (nonlinear régime), when the non- 
linearity length is much smaller than the Anderson local- 
ization length, a soliton can propagate almost unchanged 
through a disordered médium as predicted theoretically 
in 1990 H and recently verified experimentally ||. 

Thèse two important physical concepts, spatial self- 
trapping of light generated by an incohérent source in 
a homogeneous médium, and suppression of Anderson 
localisation for large-amplitude waves in spatially disor- 
dered média, both resuit from the effect of strong nonlin- 
earity. When the nonlinearity is sufRciently strong it acts 
as an effective phase-locking mechanism by producing a 
large frequency shift of the différent random-phase com- 



ponents, and thereby introducing an effective order into 
an incohérent wave packet, thus enabling the formation 
of localised structures. In other words, both phenomena 
correspond to the limit when the ratio of the nonlinearity 
length to the characteristic length of (spatial or tempo- 
ral) fluctuations is small. In the opposite limit when this 
ratio is large the wave propagation is basically linear. 

What will happen in the intermediate case when the 
length scales of nonlinearity and fluctuations become 
comparable ? It is usually believed that localised struc- 
tures would not be able to survive for such incohérent 
wave propagation and should rapidly decay. In this 
Letter we show that, at least for aperiodic inhomoge- 
neous structures, solitary waves can exist in the form 
of quasiperiodic nonlinear localised modes. As an exam- 
ple we consider second-harmonie génération (SHG) and 
nonlinear beam propagation in Fibonacci optical super- 
lattices, and demonstratc numcrically the possibility of 
spatial self-trapping of quasiperiodic waves whose enve- 
lope amplitude varies quasiperiodically, while still main- 
taining a stable, well-defined spatially localised structure, 
a quasiperiodic envelope soliton. 

We consider the interaction of a fundamcntal wave 
(FW) with the frequency u> and its second harmonie (SH) 
in a slab waveguide with quadratic (or x^) nonlinear- 
ity. Assuming the x susceptibility to be modulated 
and the nonlinearity to be of the same order as diffrac- 
tion, we write the dynamical équations in the form 
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where w(x,z) and v(x,z) are the slowly varying en- 
velopes of the FW and SH, respectively. The parame- 
ter (3 = A/c|fc w |xQ is proportional to the phase mismatch 
Afc = 2fc w — k2uj, k u and ki^ being the wave numbers 
at the two frequencies. The transverse coordinate x is 
measured in units of the input beam width xq, and the 
propagation distance z in units of the diffraction length 
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ld = %x>\kw\- The spatial modulation of the suscep- 
tibility is described by the quasi-phase-matching (QPM) 
grating function d(z). In the context of SHG, the QPM 
technique is an effective way to achieve phase matching, 
and has been studied intensively (see Réf. for a com- 
prehensive review). 

Here we consider a QPM grating produced by a 
quasiperiodic nonlinear optical superlattice. Quasiperi- 
odic optical superlattices, one-dimensional analogs of 
quasicrystals ||, are usually designed to study the ef- 
fect of Anderson localisation in the linear régime of light 
propagation. For examplc, Gcllcrmann et al. measured 
the optical transmission properties of quasiperiodic di- 
electric multilayer stacks of Si02 and Ti02 thin films and 
observed a strong suppression of the transmission Q . For 
QPM gratings, a nonlinear quasiperiodic superlattice of 
LiTaOa, in which two antiparallel ferro-electric domains 
are arranged in a Fibonacci séquence, was recently fabri- 
cated by Zhu et al. |ï^| , who measured multi-colour SHG 
with energy conversion efficiencies of ~ 5% — 20%. This 
quasiperiodic optical superlattice in LiTa03 can also bc 
used for efficient direct third harmonie génération 
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FIG. 1. (a) Slab waveguide with quasiperiodic QPM su- 
perlattice structure composed of building blocks A and B. (b) 
Numerically calculated amplitude spectrum of d(z). 

The quasiperiodic QPM gratings have two building 
blocks A and B of the length Ia and Ib, respectively, 
which are ordered in a Fibonacci séquence [Fig. |ï|(a)]. 
Each block has a domain of length ^=1 (^=1) with 
d=+l (shaded) and a domain of length Ia 2 =1(1 + rj) 
[£b 2 =Z(1 — ttj)] with d=— 1 (white). In the case of 
nonlinear QPM superlattices this corresponds to positive 
and négative ferro-electric domains, respectively. The 
spécifie détails of this type of Fibbonacci optical super- 
lattices can be found elsewhere (see, e.g., Réf. Q and 
références therein). For our simulations presented below 
we have chosen r\= 2(r — 1)/(1 + t 2 )= 0.34, where r= 
(1 + \fh)/2 is the so-called golden ratio. This means that 
the ratio of length scales is also the golden ratio, Ia/Ib= 
t. Furthermore, we have chosen Z=0.1. 

The grating function d(z), which varies between +1 
and —1 according to the Fibonacci séquence, can be ex- 
pandcd in a Fourier séries 



d(z) = d m , n e lG " 



G r , 



27t(to 



where D—tIa + Ib—0-52 for the chosen parameter values. 
Hence the spectrum is composed of sums and différences 
of the basic wavenumbers ki=2tt/D and k,2=2itt/D. 
Thèse components fill the whole Fourier space densely, 
since n\ and K2 are incommensurate. Figure ^(b) shows 
the numerically calculated Fourier spectrum G mi „. The 
lowest-order "Fibonacci modes" are clearly the most in- 
tense. From Eq. (||) and the numerically found spectrum 
we identify the six most intense modes presented in Ta- 
ble 1. The corresponding wavenumbers G rn . n are in good 
agreement with Eq. (0). 
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10 12 12 
1 1 2 3 4 
31.42 19.42 50.83 82.25 12.00 101.66 



TABLE 1. The six most intense Fibonacci modes G m ,n. 

To analyse the beam propagation and SHG in a 
quasiperiodic QPM grating one could simply average 
Eqs. (0). To lowest order this approach always yields 
a System of équations with constant mean-value coeffi- 
cients, which does not allow to describe oscillations of 
the beam amplitude and phase. However, here we wish 
to go beyond the averaged équations and consider the 
rapid large-amplitude variations of the envelope fonc- 
tions. This can be done analytically for periodic QPM 
gratings [Q. However, for the quasiperiodic gratings we 
have to résolve to numerical simulations. 

Thus we have solved Eqs. (Q) numerically with a 
second-order split-step routine, in which the linear part 
is solved with the fast-Fourier-transform (FFT) method 
and the nonlinear part, with a fourth-order Runge-Kutta 
scheme. The step-length is adapting to the local domain 
length of the QPM grating. At the input of the crystal we 
excite the fundamental beam (corresponding to unseeded 
SHG) with a Gaussian profile, 



w(x, 0) = A w e 
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v(x,0) = 0. 
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We consider the quasiperiodic QPM grating with match- 
ing to the peak at G2,3, be., /3=G2,3=82.25. First, we 
study the small-amplitude limit when a weak FW is in- 
jected with a low amplitude. Figures |^(a,b) show an ex- 
ample of the évolution of FW and SH in this effectively 
linear régime. As is clearly seem from Fig. ||(b) the SH 
wave is excited, but both beams eventually diffract. 

When the amplitude of the input beam exceeds a cer- 
tain threshold, self-focusing and localization should be 
observed for both harmonies. Figures ||(c,d) show an ex- 
ample of the évolution of a strong input FW beam, and 
its corresponding SH. Again the SH is generated, but now 
the nonlinearity is so strong that it leads to self-focusing 
and mutual self-trapping of the two fields, resulting in 
a spatially localized two-component soliton, despite the 
continuous scattering of the quasiperiodic QPM grating. 

It is important to notice that the two-component lo- 
calised beam created due to the self-trapping effect is 
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quasiperiodic by itself. As a matter of fact, after an ini- 
tial transient its amplitude oscillâtes in phase with the 
quasiperiodic QPM modulation d(z). This is illustrated 
in Fig. [| where we show in more détail the peak inten- 
sities in the asymptotic régime of the évolution. 



Fibbonacci QPM grating spectrum. Thus, low-amplitudc 
broad solitons are excited for /3-values in between peaks, 
whereas high-amplitude narrow solitons are excited when 
j3 is close to a strong peak, as shown in Fig. ||(c,d). 
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FIG. 2. (a) Diffraction of a weak FW beam with amplitude 
A w =0.25 for /3=82.25. (c) Excitation of a quasiperiodic soli- 
ton by a FW beam with amplitude A w =5 for (3—82.25. (b,d) 
Corresponding SH components. 

Since the oscillations shown in Fig. arc in phase with 
the oscillations of the QPM grating d(z), their spectra 
should be similar. This is confirmed by Fig. ^, which 
gives the spectrum of the peak intensity \w(z, 0)| 2 of the 
FW. Note that the Fibonacci peak at fc=82.25 is sup- 
pressed (or reduced) because the identical mismatch (3 
down-converts it to the dc-component. Sum and différ- 
ence wavenumbers between j3 and G m>n appear, which 
arc generated by the nonlincarity. For example, the com- 
ponent at fc=62.8 is the différence between /3=82.25 and 
G i = 19.42. 
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FIG. 3. Amplitude oscillations of the quasiperiodic soliton. 
(a),(b) Close-ups of the peak intensity \w(0, z)\ 2 of the FW 
(black) and |«(0, z)\ 2 of the SH (grey). The Fibonacci build- 
ing blocks A and B are indicated in (b) with d=l in grey 
régions, and d=— 1 in white régions. 



Our numerical results show that the quasiperiodic en- 
velope solitons can be generated for a broad range of the 
phase-mismatch (3. The amplitude and width of the soli- 
tons dépend on the effective mismatch, which is the sépa- 
ration between (3 and the nearest strong peak G m ,n in the 
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FIG. 4. Spectrum of the amplitude oscillations of the FW 
component of the quasiperiodic soliton, calculated from z~90 
to 100 in Fig. [](a). The peaks correspond to the Fibonacci 
peaks G m ,n in d(z) and sum and différence thereof with the 
mismatch /3=82.25. 



The existence of spatially localized self-trapped states 
in nonlinear quasiperiodic média should not dépend on 
the particular kind of nonlincarity. The dependence on (3 
observed here for the x^ 1 gratings is simply because the 
"real" strength of the x^ nonlinearity is inversely pro- 
portional to the phase-mismatch. In fact, it is wcll-known 
that for large values of the mismatch f3 the quadratic 
nonlinearity becomes effectively cubic |Ï3[] . Thus, our 
findings are directly applicable to nonlinear optical su- 
perlattices in cubic (or x*- 3 -*) nonlinear média. 

To analyse in more détail the transition between the 
linear (diffraction) and nonlinear (self-trapping) régimes, 
we have made a séries of careful numerical simulations. 
In Fig. H we show the transmission coefficients and the 
beam widths at the output of the crystal versus the in- 
tensity of the FW input beam, for a variety of /3-values. 
Thèse dependencies clearly illustrâtes the universality of 
the génération of localised modes for varying strength 
of nonlinearity, i.e. a quasiperiodic soliton is generated 
only for sufhcicntly high amplitudes. This is of course 
a gênerai phenomenon also observed in many nonlinear 
isotropic média. However, here the self-trapping occurs 
for quasiperiodic waves, with the quasiperiodicity being 
preserved in the variation of the amplitude of both com- 
ponents of the soliton. 

Numerical simulations for other values of the phase- 
mismatch (3 reveal the same basic property of quasiperi- 
odic self-trapping: Spatial solitons are formed in Fi- 
bonacci quadratic nonlinear slab waveguides above a 
certain power threshold, and such solitons are always 
quasiperiodic, i.e. they exhibit large-amplitude oscilla- 
tions along 2, which are composed of mixing of the two 
incommensurate Fibonacci wavenumbers K\ and «2- The 
amplitude and width of thèse solitons dépend on the dif- 



3 



ference between the phase-mismatch f3 and the nearest 
strong peak G m . n in the Fibonacci spectrum. 
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FIG. 5. (a) Fransmission of the FW, \w(0, L)/w{0, 0)| 2 , 
and the SH, \v(Q, L)/w(0, 0)| 2 vs. input peak intensity 
\w(0, 0)| 2 of the FW. (b) Output beam width vs. |to(0,0)| 2 . 
L=100, /8=Gi,2=50.83 (solid), /3=G 2 ,4=101.66 (dashed), and 
/3=G 2 , 3 =82.25 (dot-dashed). 

Finally we would like to emphasize that the phe- 
nomenon described here is qualitatively différent from 
the propagation of topological and nontopological kinks 
in disordered and quasiperiodic nonlinear média [Q . Such 
kinks can be well approximated by an effective structure- 
less particle, which either préserves identity, as in the 
case of topological kinks JÏ4|J[5| ], or decays rapidly into 
radiation |ïq| . 

In conclusion, we have analysed SHG, self-focusing, 
and nonlinear beam propagation in Fibonacci optical su- 
perlattices with a quadratic nonlinear response. We have 
predicted spatial self-trapping of quasiperiodic waves and 
the formation of quasiperiodic solitons. Such solitons 
have a localised envelope that traps the random-phase 
components through the phase and frequency locking ef- 
fect of strong nonlinearity, and whoes amplitude under- 
goes clearly détectable quasiperiodic oscillations. The re- 
sults presented here would allow to extend the concepts 
of self-localisation and self-modulation of nonlinear waves 
to a broader class of spatially inhomogeneous média, and 
can also be found in Systems of différent physical context. 
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